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Abstract 

A late time asymptotic perturbative analysis of curvature coupled com- 
plex scalar field models with accelerated cosmological expansion is carried 
out on the level of formal power series expansions. For this, algebraic ana- 
logues of the Einstein scalar field equations in Gaussian coordinates for 
space-time dimensions greater than two are postulated and formal solu- 
tions are constructed inductively and shown to be unique. The results 
obtained this way are found to be consistent with already known facts 
on the asymptotics of such models. In addition, the algebraic expansions 
are used to provide a prospect of the large time behaviour that might be 
expected of the considered models. 

1 Introduction 

A common and also remarkably consistent interpretation of current observa- 
tional data on type la supernovae distances, cosmic microwave background 
anisotropies and large scale structure surveys is that the universe is of spatially 
flat Friedmann Robertson Walker type and is currently undergoing accelerated 
expansion. However, the supposition of such a model opens a gap between the 
predicted total energy density required by a flat space, the critical density, and 
the sum of observed baryonic and dark matter density. An unknown, smooth 
component of matter with negative pressure, dubbed dark energy, is usually in- 
troduced to fill this gap of about two thirds of the critical density. For a recent 
review on dark matter and dark energy see e.g. [2B . It is noteworthy that 
there are also attempts to avert the need of dark energy at all, for instance by 
assuming Van der Waals type 0| or viscous ^U] instead of perfect matter fluids. 

Trying to fathom the nature of dark energy, myriads of models have been 
brought up for it, from the simplest, as just a positive cosmological constant [25] 
yielding an equation of state parameter w = — 1, to far richer, dynamical ones, 
where w can vary in time. Many of these models make use of various scalar 
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fields to describe dark energy, although there are not less successful approaches 
like the generalized Chaplygin gas ^B] which consider fluids that are to obey 
specific nonlinear equations of state. Models with real scalar fields minimally 
coupled to the Ricci curvature and propagating in different (non-negative) po- 
tentials, as proposed in [50], are known as quintessence and cause w lying 
in the "ordinary" range between minus one and one. If the kinetic energy of 
such a field is sufficiently small compared with the potential energy, w will be 
less than —1/3, so driving the acceleration of the universe An interest- 

ing generalization of this is spintessence where the scalar field is assumed 
to be complex or even to possess an internal O(N) symmetry ^ . The 
restrictions on quintessence can be relaxed in other ways by allowing arbitrary 
coupling to the scalar curvature [T3|, to the matter fluid pQ, or, instead of a 
potential, introducing non-standard kinetic terms 0- The latter is known as 
fc-essence. Doing so enables w to take values below —1, a case which has been 
discussed under the term phantom energy [J| . Phantom cosmologies are capable 
of exhibiting a certain type of future singularity [7J |2E] , the so called big rip [5] , 
where the scale factor diverges in finite time. 

A way of getting information about the late time asymptotics of cosmological 
models is to consider formal power series expansions of their solutions as sug- 
gested by Starobinsky in |27] and worked out in full detail by Rendall for the 
vacuum case. Recently, for example, Heinzle and Rendall calculated formal 
asymptotic expansions for quintessence in an exponential potential, whereas in 
|3] and ^31 formal series were used to investigate the big rip singularity due to 
a phantom scalar field and phantom barytropic perfect fluid respectively. The 
method presented in |21| encompasses two parts, a first algebraic one, in which 
the Einstein equations are used to obtain algebraic equations for the series' co- 
efficients that are then solved inductively, and a second analytic one, where it is 
shown that those formal series indeed approximate the cosmological solution to 
arbitrary high order for large times. Rendall |^] did also carry out the algebraic 
part for perfect fluids with linear equation of state, more specifically for con- 
stant non-negative w less than one. It is the aim of the present paper to do such 
a perturbative analysis for curvature coupled complex scalar fields. For this, 
section|2|recapitulates the Einstein scalar field equations in n + 1 decomposition 
with respect to a Gaussian coordinate system. In section[3]the formal series are 
defined as elements of the algebra of generalized finite formal Laurent series over 
a suitable ring. Section 0] then gives partial and covariant derivatives as well 
as (matrix) inverses of such series a precise meaning and postulates algebraic 
analogues of the Einstein and scalar field equations obtained before. Finally, in 
section |5l existence and uniqueness of solutions to this algebraic equations are 
shown. 

2 Motivation of the equations 

In the algebraic treatment of scalar field models in the sections following, a 
set of equations for formal series will be postulated. These equations shall be 
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motivated here as originating from a n + 1 decomposition of the dynamics in 
Gaussian coordinates. For this, let g be a Lorentzian metric on M := M x /, 
where M is a smooth manifold of dimension n > 2 and I an open non-empty 
interval within the real numbers R, such that the n + 1 decomposition has unit 
lapse and vanishing shift, i.e. 

g\U x I = g ab dx a g> da; 6 - dt®dt 

for every chart 1 ([/, x) of M, where t denotes the natural projection on I and 
g(-,to) is the Riemannian metric induced by g on the hypersurface t = to £ I. 
Then the Einstein equation G — T on M is equivalent to a system consisting of 
an evolution equation 

d t k\ = R\ + (tr k)k\ - S\ + ^—^5\ (1) 

n — 1 

as well as both an energy and a momentum constraint equation 

R-k^ + (trfc) 2 = 2p (2) 
Vik\, -V„trfc = Ja- (3) 

Here, fc a b is the scalar second fundamental form of the t-hypersurfaces in M, 
which satisfies 

d t g a b = -2fc a b, (4) 

whereas R ab and R are the Ricci tensor and scalar curvature of g. Moreover, 
the following projections of the energy momentum tensor T^ v with respect to 
the i-hypersurfaces were used |21) : 

g a b rpab. ja rpan. ^ rpnn 

Latin indices shall indicate spatial components and run from to n— 1, whereas 
Greek indices may cover the entire range, from to n. 

If (j) is a complex scalar field on M coupled to the scalar curvature R with 
strength £ G R and propagating in a smooth real potential V, then it fulfills the 
wave equation 

D(p-^R4>-24>V'(4>*4>) = (6) 
and its energy momentum tensor is given by 2 



T = i 
2 



( I - - X IV,,. ■ Y,.<> + ( 2£ - i ) ,/,,„ V„. . V. . - / V,, V„. 



+ 2ig^*W a V a ^> + ZG^pb - g»vV(4>*4>) + c.c. 



(7) 



The abbreviations □ for the covariant Laplacian A = V a V Q of M as well as 
c.c. for the complex conjugate of all the preceding terms were used. Note that 

1 Compositions with the natural projection U X I — » U are not written out explicitly. 

2 A recent and well-referenced discussion on non-minimal coupling can be found in 1111 . 
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in contrast to real scalar fields, where the potential term is usually written as 
V(cf>), it is here defined to be V((f>*4>) for having V a function of a real variable. 
The explicit occurrence of the Einstein tensor in 0) can be circumvented by 
assuming 1 — ' <f> > and introducing an effective energy momentum tensor 



_ T-ZGF<i> 

The wave equation JfjJ) is satisfied if and only if 

A</> + (tr k)d t <p - d?(/> - £R<j) - 2<j)V' 



= 



holds on M, where the ambient scalar curvature R can be written as 
R = R- 2d t tr k + hjk tj + (tr kf . 



(8) 
(9) 
(10) 



The Einstein equations are then equivalent to G — T which in turn is valid 
exactly if the (effective) evolution and constraint quantities [21] 



E\ := d t a\- 



E := atrfc- 



R\ + (tvk)a\-S\ 



i?+ (trfc) 2 + 



1 



n — 1 



-trS- 



n — 1 



C := R-kijk ij + (trfc) 2 -2p 
C a := W i k i a -V a txk-j a 



(11) 

(12) 

(13) 
(14) 



all vanish identically. In this paper, a hat shall denote the tensor's traceless 
part and a := k. The projections of the effective energy momentum tensor T 
can be written out as 
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(i - 



2£ 



(1 - 2£)V a </>* V b <j> ~ 2^* [V a V b + k\d t 
V^*V> - dt^d t <t>] + 2i(f>* [A(f> 



- (1 - 20 Wfttf + 2&* [d t V a < 



+ (tr k)d t cf> - 3U\ ~ V(4>*(/>) \S a b + c 



k a ,V l 



-ci-e0^)- x 



(15) 



(16) 



2^ - 5 ) Vi 



(trfc)a^] + y(0>) 



(17) 



4 



Finally, two useful consistency conditions can be obtained by differentiating 
the energy and momentum constraint quantities (|13l) , 1141) with respect to time 



d t C = 2(tr k)C - 2V l C t - lV, A ', + 2 ( 1 - - ) (tr k)E 



n 



-2[d tP + VJ l ~ k l3 S ij - (tr k)p] (18) 

d t C a = (tr k)C a - iv a C + ViE\ - U - V a £ 

- 9ai [d t f + VjS*' - 2^ - (tr k)f] (19) 

and identifying the matter terms in square brackets above as the temporal and 
spatial components of the covariant divergence of the effective energy momen- 
tum tensor respectively. 

3 Algebraic setup 

To investigate the asymptotics of curvature coupled scalar field models alge- 
braically, an ansatz for both the metric and the matter field in terms of formal 
power series is made. Then, the Einstein scalar field equations ©, © an d 
© provide relations on the coefficients of those series. An ansatz found to be 
suitable is of the form 



E E 



J2 « mAl (^^ (m+,s)ffi , 



mG-D s£S m iGi m , s 

where D belongs to a certain class of subsets of the reals given below and S m 
and L m ^ s are finite sets of real and natural numbers respectively, this for m £ D 
and every s € S m . Furthermore, H is a positive constant playing the role 
of an asymptotic value of the Hubble parameter, and a m ^ St i(x) is a coefficient 
independent of t. Non-vanishing terms with I ^ will be called logarithmic, 
such with s ^ oscillatory. 

In order to give the series above a precise meaning, it is now shown that they 
can be taken as elements of a specific algebra providing the structure required. 
For this, let R be a complex algebra 3 and R{(X)) the set of functions from 
M. to R whose support is discrete 4 and bounded from below. Then this set 
together with point-wise addition, point-wise exterior multiplication and the 
Cauchy product forms the complex algebra of generalized finite formal Laurent 
series in one variable X over R, where the supposition on the support of the 
functions guarantees that the sums in the inner multiplication are all finite and 
therefore well-defined. Note that for complex k £ C and f,g £ R((X)) the 
support of nf is a subset of the support of /, that of / + g is a subset of the 
union of the supports of / and g, whereas the support of the product fg lies 
within the point-wise sum of the supports of / and g. As usual, X a denotes the 

3 Throughout this paper, an algebra is considered associative, commutative and unitary. 
4 A subset A of a topological space X is called discrete iff A has no limit points in X. 
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map for which X a (m) = 1r, the identity in R, for m — a and zero otherwise. 
By virtue of the identity 1r{{x)) = X° in R((X)) the algebra R is identified with 
the subalgebra R1 R ^ X )) m In addition, a function taking values in R 

defined on a subset of M is always identified with its trivial extension, where 
zero is assigned to points outside the function's domain. It will be useful to say 
a / S R((X}} being of order not less than a, written / = 0(X a ), iff a < supp /, 
being of order greater than a, written / = o(A Q ), iff a < supp/, i.e. non-zero 
coefficients occur not before or even after a € M respectively. 

There are three subalgebras of R((X}) which will be used further on. Firstly, 
the polynomials R[X] over R encompassing exactly those series supported on a 
finite subset in the natural numbers, secondly, the power series over R, 

obtained by allowing series supported on any set of naturals, and thirdly, the 
generalized polynomials R(X) over R, whose support is a finite subset in the 
real numbers. With this, and by noting that the set of smooth complex valued 
functions C°° (M, C) on the n-manifold M possesses the structure of a complex 
algebra, it is possible to define an algebra C := C°° (M, C) [Z] (Y) ((X)) over C 
in three variables X, Y and Z successively, i.e. C = ((C°°(M, C)[Z])(Y))((X)). 
The identification of these variables with distinct elements of C follows from 
the outline in the previous paragraph. An element of the algebra C can be 
regarded as a series of the form © when the variable X plays the role of the 
exponential factor e~ Ht , Y that of the oscillatory factor e~ lHt and Z that of 
the logarithmic factor t. The coefficients are smooth complex functions on M. 
It is worth remarking that this definition of the formal series is totally intrinsic 
to M, the appearance of t in Q is merely symbolic and intended to suggest 
arithmetic rules rather than being a reference to the ambient manifold M. To 
have the setup complete, let S u denote the subalgebra of exactly those elements 
of £, which are of order not less than zero and have neither oscillatory nor 
logarithmic terms occurring before Using the above identifications, the 

absence of those terms can also be expressed as the restriction of the series to 
the open interval ]-oo, u[ lying in C°°(M, C)((X)). 

4 Algebraic Einstein scalar field equations 

After these preparations, it is possible to postulate Einstein and scalar field 
equations in terms of formal series. For a / € C, formal spatial derivatives 
shall be defined to act point-wise, (9 a /)m,s,i : = da(f)m,s,h an d the formal time 
derivative according to (<9 t /) m , s ,j := -(m + is)H(f) miS j + (I + l)(/) m ,s,J+i> 
77i, s, Z G R. This imitates the effect the partial derivatives would have when 
applied term-wise to an ansatz of the form @. Here, the notation (f) m ,s,i for 
/(m)(s)(Z) € C°°(M, C) is used. Those formal derivatives are derivations on the 
algebra £ commuting pairwise and leaving the subalgebras 5„,u£l, invariant. 
To have a notion of reality for formal series, define complex conjugation as a 
ring homomorphism / h-> /* by (f*) m , s ,i ■= U)* m .- S ,i for a11 m,s,l eR and say 
that / is real iff /* = /. 

Fix a chart on M which, for the sake of simplicity, may be global and let 
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{9ab) be a real symmetric family in S U X~ V for some u > and v G M with 
(g a b)-v,o,o (the components of) a Ricmannian metric. It can be shown that there 
exists exactly one real symmetric family (g ab ) in S U X V with g a ig lb = 5 a b and 
g a% gib = S a b . Thus, together with the formal partial derivatives above, algebraic 
equivalents for the components of the Levi-Civita connection, the Ricci tensor 
and hence for covariant derivatives and scalar curvature of M can be obtained 
by setting 

= ^g lc (d a g lb + d b g a i ~ dig a b) and 
Rab = d^lb-daT^ + K^-Tl^. 

In addition, equation is used to define k a b algebraically and, in this context, 
also intrinsically to M. 

Having the formal series for the geometric quantities at hand, algebraic rep- 
resentations of the components of the energy momentum tensor can be found. 
For this, let £, //, A £ R be arbitrary values for the coupling constant, field mass 
and cosmological constant respectively and assume H > 0. A formal scalar field 
tf> £ jC of order higher than zero and a smooth self interaction potential shall 
be given such that V(0) = A and V'(0) = /i 2 /2, i.e. the first two coefficients 
of the McLaurin expansion may be denominated explicitly A formal equivalent 
of V(cj)*(j)) is naturally obtained by applying the substitution homomorphism 
induced by (j)*^ to the McLaurin expansion of V, hence setting 

w ■■= E -,(d n v)(o)(4>*<t>r 

*■ — ' n\ 

n=0 

which is well-defined because of (jfcj) = o(l) and so the sum is point-wise finite. 
Substitution into the expansion of V shall yield W' in the same manner. Also 
due to (f>*(f> = o(l), the element 1 — £</>*</> has a multiplicative inverse in C which 
is exactly of order zero. Now the components of the effective energy momentum 
tensor l|15|) . (|lfc>|> and (|17ll . the evolution and constraint quantities (flip . (|12|) . 
(|13|) and (|14|l as well as the scalar field equation make sense in terms of 
formal series and therefore lead to the following definition^ It is convenient to 
introduce the abbreviations 

[n 2 1 'A - ^n(n + 1)] H 2 

f n/2 for [i 2 > K 

\ n/2 - \]K — [i 2 /H for ^ < K 

f n/2 for n 2 > K 

\ n/2 + y/K - t i 2 /H for [i 2 <K 

here. 

Definition 1. Let (g a b) be a real symmetric family in S U X~ V for some u > 
0. v G BL with (g a b)-v,o,o a Riemannian metric and <f> G C of order greater 
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than zero. Then the pair (gabj'P) * s sa id t° be a solution of the (algebraic) 
Einstein or scalar field equations at relative order m £ M. iff the evolution and 
constraint quantities vanish at m, so (E a b , E,C,C a )(m) = 0, or holds at 
m+k\ respectively. It is called a solution of the Einstein or scalar field equations 
iff it is a solution of the corresponding equation at every relative order m £ K. 

Furthermore, the concept of relative orders shall be extended to the metric 
and field by defining the coefficients of g ab and 4> a t relative order m to be 
gab( m — 2) and <p(m + k\) respectively. 

Simple but useful necessary conditions for the existence of algebraic solu- 
tions as well as relations of the asymptotic Hubble constant H or the principal 
frequency of field oscillations to the cosmological constant A or field mass fi, 
respectively, can be obtained by focusing attention to the lowest order terms in 
the series for both metric and field. Let (g a b, </>) be a solution of the Einstein 
scalar field equations. Prescribing the metric g a b to lie within S n X~ 2 and have 
coefficients (<? a b)-2,o,o that make up a Riemannian metric, implies, as mentioned 
above, g ab £ S n X 2 which in turn requires the components of the connection r^ b 
and Ricci tensor R ab being real and in S n . This in turn leads to R a b , R £ S n X 2 
both real, whereas equation Q identifies k ab £ S n X~ 2 real and so either k a b 
and tr k as real elements of S n . The same equation (0J also determines their 
lowest order coefficients entirely, 

(k a b ) , ,o = -H6 a b ; (<rY)o,o,o = 0; (tr fc)o, , = -nH. (20) 

Further, assume the field (j) being in Sk 2 and not equal to 0, then the mini- 
mum r of the support of <p exists and the scalar field equation at r reads 

(I + l)(i + 2)(£) r , <M+2 + (I + l)(n - 2r - 2is)H(<t>) r>s ,i + i 

+ {n 2 - s 2 H 2 + [r(r - n) + £77(77 + 1)] H 2 + is(2r - n)H 2 } (<p) r , Sl i = 



for all s, I £ K. From this, it can be seen that <j> is actually an element of 

&2 — ' 



Sh-j-k, X kl and that the assertion 



0(r) £ C°°(M,C) for /j 2 < K 
4>{r) £ C°° (M, C) + C°° (M, C)Z for /j 2 = K 
4>{r) £ C°° (M, C)Y-" + C°° (M, C)Y" for /j 2 > K 

holds, i.e. for subcritical field masses there are neither logarithmic nor oscil- 
latory initial terms, for critical field masses at most one logarithmic but no 
oscillatory initial terms and for supercritical field masses no logarithmic but in 
general two oscillatory initial terms at frequency u>H present, were 



._ / VM 2 - K/H for 11 2 >K . . 

) for fi 2 <K. [Zi) 

Using the lowest order coefficients for k a b known from J5DJ, the scalar field equa- 
tion (0 can be written in a form suitable for calculating the field inductively, 
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namely for every positive relative order m > in case of subcritical or critical 
field masses fj, 2 < K as 

(I + + 2)(<f>) m+kuSt i +2 -(1 + l)(2m + fei - fc 2 + 2is)i3"(^) m+fel>s , i+ i 

+ (m + is)(m + fei - fc 2 + is)H 2 (4>) m+kltS> i = LOT m+kl , s ,i, (22) 

in case of supercritical field masses /i 2 > K as 

(i + + 2)(^) m+fel , s>i+ 2 - 2(i + l)(m + is)H(4>) m+kuS ,i+i 

+ (m + is- iu)(m + is + iuj)H 2 (4>) m+kl ^j = LOT m+feliSj /, (23) 

where in LOT m -|_fc liS) j only field coefficients of relative order lower than m and 
metric coefficients of relative order lower than or equal to m occur. 

From the fact that <fi is in Sk 2 - kl X kl it follows that <p*<P, W and W 1 lie within 
S n , so do the components of the effective energy momentum tensor S a b ,j a ,p 
and with them, finally, all the evolution and constraint quantities E a b , E, C, C a . 
Furthermore, according to the equations i|f 5|l . (|f 6fl and (|f 7|l . 

(S Q Jo,o,o = ~ AS \ ; (J'a)o,o,o = 0; (p)o,o,o = A (24) 

is implied, and hence the validity of the Einstein equations to relative order zero 
requires in particular A > and 



y n(n — 1) 

Their validity up to but not including fco := min{2,2fci} yields the vanishing 
of g a b and hence k a b on the open interval ]— 2, ko — 2[ , of g ab on the interval 
]2, fco + 2[ , thus of k a b , T c ab and R ab on the interval ]0, fco[ an d so of R a b on 
]2, k + 2[ by induction. 

Taking the last statement above into account, it is possible to express the 
evolution and constraint quantities in a coefficient-wise manner for relative or- 
ders m > fco as 

(£ a b )m, s ,i = (n-m- is)H(a a b ) m>St i + {I + l)(a a b ) m , s ,l+i - {R a b ) m ,s,i 

- E E E (<7 a b)p,uAtek) q , v ,x + (S a b ) m , sJ (26) 

p+q=m u+v=s k+X=1 

{E) m ,s,i = (2n - m - zs)_ff(tr fc) miSi ; + (I + l)(tr fc) miS , i+ i - (R) m .s,i 

~ E E E ( trk )p.uA tTk )q,v,\ 
p+q=m u+v=s k+X=1 

- [(tr S) m ,.,i + n{p) m ,,,i] /(n - 1) (27) 
(COm.s.J = ~(2n - 2)fr(trfc) m , s ,i + (R) m ,s,i ~ 2(p)m,s,l 

+ E E E 

p+q=m u+v=s k+X=1 

+ (trfc) p ,„, re (trfc)g )W>A ] (28) 
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(C a )m,s,l — ^ i(^'o)m,s,l _ V o^r^m^,! - (i a )m,s,l 

- E E E [(riiW(^)*M 

p+g=m u+-u=s K+A=i 

-(rkW^w], (29) 

where the indices p and q are never less than ko and V° denotes the Levi-Civita 
connection induced on M by the Riemannian metric (<7ab)-2.o.o- For the same 
relative orders, equation Q reads 

(m + is)H(g ab ) m ^ 2 .s,i - (I + l)(g a b)m-2,s,l+l = 

2(9ai)-2fi,o{k\) m ,sJ + 2 E E E i9ai)p-2,uA k \)q,v,\- (30) 
p+q=m u+u=s K +\=l 

Assuming validity of the Einstein scalar field equations up to but not in- 
cluding relative order m > k , the consistency conditions 118fl and (|19|) simplify 
to 

(2n- m - is)H{C) m}St i + (I + l)(C) m ,,, l+1 = -2(n - l)H(E) m ,,,i (31) 

and 

(n - m - is)H(C a ) m ,s,i + (I + l)(C a ) m , s ,i+i = 

-\^ a {C) m , s>l + V\(E\) m . Sil - (l - V„(E) m , s ,,. (32) 

5 Construction of the solution 

Using the relations for the coefficients of metric and scalar field derived in the 
previous section, existence and uniqueness of algebraic solutions of the Einstein 
scalar field equations can be shown by (transfinite) induction. After making 
some considerations about the support of the required formal series, initial so- 
lutions are constructed up to relative order n, exclusive, which by themselves 
will allow stating a necessary and sufficient condition for the existence of solu- 
tions up to beyond relative order n. For this, let a coupling constant £ S R, a 
cosmological constant A > 0, a field mass fi € M and a smooth real valued self 
interaction potential V, defined on an interval containing zero, with V(0) = A, 
V'(0) = £t 2 /2 be given, define H and u> by l|25|l and l|21|l respectively and assume 
fj, 2 > — £n(n + 1)H 2 . The condition on the field mass ensures that k\ is positive 
and therefore the constructed field will be of order greater than zero. 

Let D be the smallest +-stable discrete non-negative set containing 0, 1, k\ 
and k2 — k\, where +- stable means the set is invariant under point-wise addition, 
i.e. D + D = D. For arbitrary g a b and </> both in C suppose (g a b)-2,o,o to be a 
Riemannian metric, the support of g a b to be contained in D — 2 and that of </> to 
be contained in D+k%. Then, the support of k a b is a subset of D — 2, the support 
of g ab one of D + 2 and so the supports of all the terms in the components of the 
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effective energy momentum tensor (JSJ , in the evolution and constraint quantities 
Qllfl. H12p. (|13() . I|14l) and hence of those quantities themselves lie within D. The 
same is true for the left hand side of the scalar field equation This means 
that it suffices to verify the validity of the Einstein scalar field equations on the 
set D. Since D is well-ordered by the induced order of IR, transfinite induction 
is available immediately. Of course, because D is at most countable, ordinary 
induction could also be used after identifying D with its isomorphic ordinal. 

As mentioned above, existence and uniqueness of initial solutions will now 
be proved. More precisely: 

Lemma 2. Let A ab £ C°°(M) be the components of a Riemannian metric and 
4>o,(j)i £ C°°(M, C) smooth complex valued functions on M. Then there exists a 
solution (g ab , 4 1 ) of the Einstein scalar field equations up to relative order n > 2, 
exclusive, with g ab & S n X~ 2 , <fi £ Sk 2 -k 1 X kl such that 

(Gl) g ab (-2) = A ab 

{00 for 
<Po + <PiZ for 
<j) Y^ + <f>xY w for 

(P2) <p(k 2 ) £ 0! + C°°(M, C)Z for ^ 2 < K 
are valid. The solution is unique up to, but not including, relative order n. 

Proof. The construction is done by induction over relative orders in the set D n 
]— oo, n[ using the following argument: For a < m < n let g a b ■ ]— oo, m — 2[ — > 
C°°(M) and 4> : ]-oo,m+fci[ C°°(M,C)[Z](Y) be constructed such that 
(gab) is symmetric, the Einstein scalar field equations are satisfied up to relative 
order m, exclusive, and that for positive m the conditions (Gl), (PI), for m + 
ki > fc 2 additionally (P2), hold. Moreover, the support of g ab may be finite, 
having exactly the point —2 lying below k — 2, the support of (f> be finite as 
well, having no points lying below k\. Finally, for all t less than both m + ki and 
/c2> <t>{t) may belong to the subalgebra C°°(M, C). Then the maps are extended 
to relative order m dependent on the value of m: 

m = : Set g a b{-2) := A ab , <p(ki) := 0o for ^ 2 > K , H k i) ■= 4>o + 4>iZ for 
= K and 0(fci) := 4> Y- U + faY" for [i 2 > K, then (Gl) and (PI) are 
satisfied, the Einstein equations hold at relative order due to the choice 
of H, the scalar field equation holds there due to the choices of u> and k\. 

< m < fco : Setting g a b(jn — 2) := fulfills the Einstein equations at relative 
order m because fco is less than or equal to both 2ki and 2. If /i 2 > K 
or fi 2 < K and m + fci 7^ k 2 , the relations 12 -ill or (|22|l yield a 0(m + 
fci ) in C°° (M, C) [Z] (Y) such that the scalar field equation is satisfied at 
relative order m. (The occurring linear systems can, for example, be solved 
successively by starting at sufficiently large values of I and descending 
towards zero.) For m + fci < k 2 4>(m + ki) actually lies within C°°(M, C). 
If, on the other hand, /j 2 < K and m + fci = k 2 , then relation (|22|) implies 



H 2 <K 
H 2 =K 
fi 2 >K 
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the existence of a \ € C°°(M, C) such that the scalar held equation is 
solved at relative order m by setting </>(m + fcx) := 4>o + X-^- This, in 
particular, suffices (P2). 

k <m <n : Using equations l(2^)l. (|27|) and l|3Ujl. g a6 (m - 2) € C°°(M) are 
found such that (g a b) is symmetric and cancels the evolution quantities 
E a b (m) and E(m) at relative order m. Now the consistency condition l|31|l 
shows that the Hamiltonian constraint quantity C{m) at relative order 
m vanishes and 1-iLl subsequently implies the same for the momentum 
constraint quantities C a (m). Thus, the Einstein equations are solved at 
relative order m. The extension of 4> can be carried out exactly as in the 
previous case: If /i 2 > K or < K and m + k\ ^ A)2, the relations l|23[l 
or yield a cf>(m + hi) in C°° (M, C) [Z] {Y) such that the scalar held 
equation is satisfied at relative order m. For m+ki < ki <f>(m+ki) actually 
lies within C°° (M, C) . If, on the other hand, /i 2 < K and m+fci = /c2 , then 
relation ll'L'l) implies the existence of a % S C°° (M, C) such that the scalar 
held equation is solved at relative order m by setting (f>(m+ki) := cf>o+xZ . 
This, in particular, suffices (P2). 

Uniqueness up to but not including relative order n follows with the same argu- 
ment by taking advantage of the necessary conditions derived in section 0] □ 

The condition for the existence of a solution of the algebraic Einstein scalar 
held equations is stated by means of a scalar £ and a 1-form Z b . For defining 
those, assume A ab being a Riemannian metric and 0o> 4>i smooth complex valued 
functions on M . By virtue of Lemma|21and equations (|26|l . J57J), l|30|) there exists 
a solution {g a b, 4>) of the Einstein scalar held equations up to relative order n, 
exclusive, with (Gl), (PI), (P2) valid and such that the coefficients (E a b ) n fij 
and (E) nj o,i+i of the evolution quantities vanish too for all natural I. It can 
then be shown that the right hand sides of 

C(A,fo,&) := ^( 5!J ) n _ 2A o-^(B)„,o,o (33) 

Z b (A, <fo, cf>i) := V% jb )„- 2 ,o,o - ~4th (V°b(-E)n,o,o + nH(C b ) n , ,o) (34) 

are independent of the choice of such a solution (g a b, <f>) and therefore the 
functions £ and Z b are well-dehned. Now the necessary part of the condi- 
tion is obvious: If (g a b,<t>) is a solution of the Einstein scalar held equations 
for which (Gl), (PI) and (P2) holds, then A«(9ij)»-2,o,o = C(A</>o,<M and 
V° (<?ib)n-2,o,o = Zb(A, 4>q, 4>i)- The sufficient part is contained in the next 
theorem. 

Theorem 3. Let A ab £ C°°(M) be the components of a Riemannian metric, 
B a i, £ C°°(M) those of a symmetric tensor and (f>o,<t>i € C°°(M, C) smooth 
complex valued functions on M with 

A l] B l3 = £(A,<h,<k) and V° l B lb = Z b (A,0 o ,0i). 
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Then there exists exactly one solution (g ab ,4>) of the Einstein scalar field equa- 
tions with g ab £ S n X~ 2 , £ Sk 2 -kiX kl and such that 



(Gl) 


(<?as6)-2,0,0 


= A ab 








(G2) 


(<?ab)n-2,0,0 


= B ab 












( 00 


for 


I' 2 


< K 


(PI) 


<f>(kl) 


= I 0o + <f>iZ 


for 


l< 2 


= K 






{ 4> Y-" + faY" 


for 


M 2 


> K 



(P2) 0(fc 2 ) e 0i +C°°(M, C)Z for < K 



are valid altogether. 

Proof. The assertion is proved by induction over relative orders in D. By virtue 
of Lemma it suffices to consider relative orders greater than or equal to n. So, 
for a m > n let g ab : ]-oo,m- 2[ -> C°°(M)[Z](F) and : ]-oo,m+ k x \ -> 
C°°(M,C)[Z](y) be constructed such that (</ a b) is real symmetric, the Einstein 
scalar field equations are satisfied up to relative order m, exclusive, and the 
conditions (Gl), (PI), (P2), and for m > n also (G2) hold. Moreover, the 
support of g ab may be finite, having exactly the point —2 lying below k — 2, 
the support of be finite as well, having no points lying below k\. Finally, for 
all t less than n — 2, g a b{t) may belong to C°°(M), as for t < k 2 , <j){t) may be 
within the subalgebra C°°(M, C). Then the metric is extended to relative order 
m as follows: 

m = n : The equations 126|) , (|27|l together with (|30|) ensure the existence of real 
g a b{n- 2) £ C°° (M, C) [Z] (Y > with {g a b)n-2,o,o = B ab such that (g ab ) is 
symmetric and causes E a b = and E(n) £ C°°(M, C), i.e. all but the 
coefficient (-E) n ,o,o necessarily vanish. Due to the trace condition implied 
on B, this coefficient turns out to be zero too, so E(n) = and by consis- 
tency equation l|31l) C(n) = is obtained. In the same manner, relation 
(13211 firstly implies C a (n) £ C°° (M, C) and the divergence condition on B 
then guarantees that C a (n) actually vanishes. So, the Einstein equations 
are satisfied at relative order n. 

m = 2n : In this case, the relations <|26|) . (|28|l and (|30|l give a real symmetric 
family of g ab {m-2) £ C°°(M, C)[Z](Y) so that E a b (m) = and C(m) = 
are fulfilled. The consistency conditions l|31|) and (|32|) then show E(m) = 
and C a (m) = successively and so the validity of the Einstein equations 
at relative order m. 

m ^ {n, 2n} : Here, the relations (|26ll . I|27|) and 13U|) yield a real symmetric 
family of g ab (m-2) £ C°° (M, C) [Z] (Y) such that E a b {m) =0 and E(m) = 
are satisfied. Once again, the consistency conditions l|3*D) and (|3*2")l are 
used to show C (m) = and C a (m) — successively and so the validity of 
the Einstein equations at relative order m. 

The continuation of the scalar field can, in any case, be done using the relations 
<|2"2"|) and respectively, to provide a coefficient 0(m+fci) £ C°°(M, C)[Z](Y) 
such that the scalar field equation holds at relative order m too. 
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The uniqueness of such a solution follows with the same inductive argument 
from the presuppositions of the theorem by making use of the uniqueness of the 
initial solution as stated in Lemma [3 □ 

It is worth noting that for non-subcritical held masses the supports of the 
series become especially simple because then k\ = fea = n/2 and so the set D 
consists of natural or halves of natural numbers only, depending on whether n 
is even or odd, respectively. Oscillatory terms by their part occur exclusively 
at integer multiples of w. If the field mass is not supercritical, there are no 
oscillatory terms at all. 

The restriction to masses /x 2 > — + 1)H 2 and thus to fields of order 
greater than zero is stringent for the treatment outlined in this section since just 
by allowing equality the right hand side of l|25|l becomes, in general, dependent 
on the field and so the Einstein equations cannot necessarily be satisfied to 
relative order zero with a constant H. 

6 Discussion 

In this paper, existence and uniqueness of formal power series solutions for com- 
plex scalar field models with arbitrary coupling to the scalar curvature is proven 
for a large class of self interaction potentials. The class includes potentials which, 
written as a function of the absolute field amplitude, possess a positive mini- 
mum at zero. The algebraic analogues of the Einstein scalar field equations 
were obtained by considering the model in n + 1 decomposition with respect 
to Gaussian coordinates and postulating those relations to hold on a suitable 
algebra of formal series endowed with all the necessary operations. This yielded 
algebraic equations for the series' coefficients which could be solved inductively. 
The series used here are generalizations of those introduced by Rcndall 21 for 
a similar treatment of the vacuum case and perfect fluid models since, apart 
from oscillatory terms, both logarithmic terms as well as non-integer powers of 
the exponential factor e~ Ht may occur. 

The lower mass limit that is imposed in the case of non-positive coupling 
constants can be interpreted as a condition for the field to decay exponentially 
in the future. For a more specific classification of the asymptotics, a critical field 
mass \f~K = ^/?i 2 /4 — £n(n + 1)H is found that distinguishes solutions based 
on the occurrence of logarithmic or oscillatory terms in the leading order k\ 
of the field's expansion. For subcritical masses, /i 2 < K, neither logarithmic 
nor oscillatory terms are present in leading order, whereas for critical masses 
/i 2 = K , the lowest logarithmic term can appear already initially. Both corre- 
spond to fields that would fall off exponentially with increasing time and this 
asymptotically at least like e~ klHt in the former, at least like e~" Ht for all /3 less 
than ki in the latter case. Thereby, the exponent k% increases from arbitrary 
small positive values for the lowest permitted masses to n/2 for the critical field 
mass. Supercritical masses on the other hand allow oscillatory terms in the 
series to show up, at frequency ±luH — ±yj u 2 — K in leading order k\ = n/2 
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and together with overtones thereof in higher orders. This has an interpretation 
as a field that could oscillate asymptotically at frequency ojH and whose am- 
plitude would decay at least like e - nHt / 2 for large times. The higher frequency 
oscillations would be attenuated exponentially with respect to the amplitude of 
the first harmonic. 

It should be emphasized that the connection between certain coefficients 
in the formal series and the asymptotics of a scalar field as made above is 
just a prospect of the large time behaviour the solutions might be expected 
to show. There is no assertion that the series would in fact be asymptotic 
in the sense of theorem 5 of or even convergent, when t is considered a 
time coordinate as in section [21 Drawing this analogy nonetheless, the results 
obtained here are found to coincide with the asymptotics of scalar field models 
derived analytically by Rendall 22 1 2.3]. In particular, the expression for the late 
time limit of the Hubble parameter (|25|) or the asymptotic value —nH for the 
trace of the second fundamental form are identical. Additionally, the condition 
on the second derivative of the potential necessary for an oscillating decay of 
the field agrees with fi 2 > K . It might also be interesting that recent SN la data 
suggests, although not by itself significantly enough to rule out even ACDM, 
the existence of oscillations in the cosmological scale factor |171 118|. In |17|. 
Lazkoz et al. achieved their best fit with a purely empirical oscillatory model. 
The observed decay of the oscillation's amplitude with the third power of the 
scale factor is compatible with the first possible occurrence of oscillatory terms 
in the expansion of the metric found here. 

Part of this work was done in the course of a diploma thesis written at the 
University of Berne, Switzerland, supervised by Prof. Dr. P. Hajfcek, in 2004. 
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